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Тригонометрические вычисления 
в системах компьютерной математики 
учебного назначения 


В работе изложены алгоритмы тригонометрических вычислений, основанные на построениях канонических 
форм. Эти алгоритмы используются при программировании стандартных задач школьного курса тригоно- 
метрии: задачи упрощения выражения, доказательства тождества, решения уравнения и неравенства. 


Введение 


Разработка алгоритмов алгебраических вычислений — одна из основных задач, 
возникающих при разработке систем компьютерной математики. Математическая 
модель задачи построения алгоритмов алгебраических вычислений — многосортные 
алгебраические системы (МАС). Практика разработки даже достаточно простых систем 
компьютерной математики учебного назначения (СКМУН) [1-3] показала, что реали- 
зация алгебраических вычислений нуждается в тщательном предварительном проекти- 
ровании МАС путем разработки иерархий сортов МАС и спецификаций интерпрета- 
торов многосортных алгебраических операций [4], [5]. 

Тригонометрия как учебная дисциплина занимает особое место. Тригонометрию 
изучают как отдельный раздел школьной алгебры и широко используют в приложениях. 
Поэтому и СКМУН должны поддерживать тригонометрические вычисления. Для реали- 
зации вычислений мы используем систему алгебраического программирования АР [6-9]. 

Целью данной работы является разработка методов тригонометрических вычисле- 
ний, основанных на построении канонических форм тригонометрических выражений. 


1 Каноническая форма целых тригонометрических 
выражений 


Определение 1. Зафиксируем множество переменных Гата Ме и будем считать, 
что на ГайаЫМе введен линейный порядок. Целыми тригонометрическими атомами 
(функциями) будем называть выражения вида 5т(Х), Со5(У), где Х,У — линейные ком- 
бинации переменных и константы РТ с рациональными коэффициентами (аргументы 
тригонометрических функций). Векторное пространство аргументов тригонометри- 
ческих атомов обозначим через ТАге. Векторное пространство линейных комби- 
наций переменных с рациональными коэффициентами обозначим через ТАте0. 

Тригонометрическим мономом (Т-мономом) будем называть произведение несколь- 
ких тригонометрических атомов. Полугруппу Т-мономов обозначим через ТМоп. 
Целым тригонометрическим выражением (Т-полиномом) будем называть линейную ком- 
бинацию Т-мономов с коэффициентами из некоторого расширения ТСое/ поля рацио- 
нальных чисел Ка! . На множестве Т-полиномов определены операции кольца и умно- 
жения на константу, то есть операции линейной алгебры. Эту алгебру обозначим 
через ТРо/!. 
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Замечание. Точное определение Т-полинома в терминах алгебраического про- 
граммирования задается спецификациями алгебр аргументов, тригонометрических 
функций, поля коэффициентов и собственно алгебры Т-полиномов. 

Конструктивная реализация алгебры Т-полиномов. Рассмотрим атом ЕКа, 
Ь, Х), который определяется формулой 


ЕКа,Ь,Х) =аСо5(Х)-+Ь5т(Х),а,Б Е ТСоеГ, Х Е ТАгэ0. (1) 
Для Х е САге пусть 
Х=есх+...+с.х, Сл, с, Е Ваь х, е ГТапаЫе, х>х,>...>х,. (2) 


Введем такие обозначения: Геа@СГ(Х) = с,,ЕгееС = с. 

Теорема 1. Любой Т-полином Р(Х)над полем ТСоеГ можно представить в 
виде суммы Е1-атомов 

РХ)=ЕКкКа,6,Х)+..+ЕКа„,б,Х,), (3) 

где а‚,Б, принадлежат некоторому алгебраическому расширению поля ТСое], 
причем, если ГеаАСУ(Х ,) > 0, Егее СГ(Х ‚)=0, такое представление единственно с 
точностью до перестановки Е 1-атомов РКа „6 »Х,), ЛЕЪ..т. 

Доказательство основано на формуле понижения степени 


ас- = ты ет Ве -аа 


ЕКа,Ь,Х). ЕКс,а,У) = ЕК ++ ВХ -У). 


Легко видеть, что (3) — сокращенная форма записи полинома Фурье 


Х,)=У`а,Со5(Х ,)+Ь,5т(Х,). 


= 


Ека,в,Х,)+...+ЕКа„,Ь 


Таким образом, 
зт(х) = РИКО, х), со$(х) = ЕКЬЬ0,х). (4) 

Расширение поля ТСое{ заключается в присоединении алгебраических чисел 
Со5(кл/п),5т(кл/п) к полю КаЁ. Вычисления в этих расширениях будут рассмот- 
рены ниже. Множество выражений вида (3) обозначим Е1Ро/. Алгебра Е1Ро[ изо- 
морфна алгебре ТРо/. 

Основная идея тригонометрических вычислений в алгебре Т-полиномов состоит в 
изоморфном переходе от выражения в сигнатуре ТРо/ к выражению в сигнатуре 
ЕРо[ (4), построению канонической формы этого выражения в Е1Ро/ и обратному 
переходу (1) в сигнатуру ТРо/. 


2 Каноническая форма рациональных 
тригонометрических выражений 


Определение 2. Пусть Ё,С - Т-полиномы от одной переменной х. Рациональ- 
Е 
ное выражение Н = а будем называть рациональным тригонометрическим выражением 


(ТВ-выражением). 
Приведя Ё,С к канонической формев Е1Ро/, получим: 


ЕКа,,Б,,Кх)+...+ЕКа,В,х)+ а 


Н(х)= . 
ЕКс,,а»„тх)+...+ ЕК, Ах) + су 


т? 
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Легко показать, что это представление не является каноническим. Зафиксируем, во- 
первых, поле коэффициентов ТСое/ поля ТК-выражений ТКаЁ. Во-вторых, определим 


и зафиксируем структуру этого расширения над полем ТСое/. Пусть ТСое{ = О- поле ра- 
циональных чисел или некоторое его вещественное алгебраическое расширение (напри- 
мер, поле Кайа= ВаТ(/2,\3,,..„|р..). Тогда поле ТКа!=О(ЕКа,Ь,Кх)),а,6 Е ТСоерК Е М. 
Рассмотрим последовательность расширений 


ТСоеГг с ТСоеГ (Соз (х)) < ТСоер (Соз (х))($т (х)) , (5) 


присоединяя к Тсое/ последовательно элементы Со5(х), 5т(х). Для краткости введем 
обозначение ТСоей =Ё,,Ё(соз(х))=Е,Е (1(х)) = Е, = ТКаЕ. Тогда расширение Ё, < Ё, 
является трансцендентным, а расширение А! < Ё, — алгебраическим. Поэтому элемент 


Н поля Ё, можно представить в виде Н =. где /,2еН[5Их)|. Многочлены Х, & 
2(х 


записаны в виде суммы мономов от степеней соз’(х), формулами понижения степени 
могут быть переписаны в многочлены от Со5( Хх). 


(>) = а Со5(Ёх) +...+а, Со5(х)+а,, &(х) =ВСо5(тх) +...+Ь, Со5(х)-+Ь 


т > 


поэтому элемент Н ЕЁ, можно представить в виде: 
Н- 1) ы аСо5(х) +...+а, (Со5(х)+а, 
2(х) ВБСо5(тх)-+...+Ь,Со5(х)+Ь, 


т-—1 
При этом Т-полиномы /,2 взаимно просты: 2с4(/,2)еЁ,. Это представление 


(6) 


однозначно с точностью до коэффициентов: для каноничности знаменатель дроби (6) 
нужно нормализовать, разделив числитель и знаменатель на [с(®) =Вь: 


= а Со5(кх) +...+а, (Со5(х)+а, — _ “4, в _ в 1) 
2(х)  Со(то+..+Ь Со5(х+Ь, ” " В, —. В, | 


2 2. 
Полученное представление уже канонично. Далее, поскольку Со5’(х)+ 5 т” (х) =1, 


степень расширения РА, < Р, равна 2. Поэтому элемент поля ЁР› можно представить 
в виде линейного двучлена: 


НЕЕ > Н=С,5(х)+С,С,,С.ЕЕ. (8) 


Итак, справедлива. 
Теорема 2. Канонической формой ТВ-выражения является его представление в виде 


Н = О аа ло) › Л›8› Л» 81 е Е = Е (соз(х)), 
8 (х) & (>) 


8с4(4, 8) =1, вса( Л, 81) =1, (8) =1, (с(в1) =1 


Для решения задачи осталось привести алгоритм обращения целого тригоно- 
метрического выражения, то есть алгоритм построения канонической формы выраже- 
ния вида 


(9) 


‚ Го, <) ЕЕ 


ЕЕ ЕВЕ 
9 тер б9+ #69 а) 
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Положим А= /(х), В = е(х), и = Со5(х), у = 5т(х). Тогда 


=УР+О, РОЕН. (11) 


Уравнение (11) относительно Р.О решим в поле Е; =Ё\(и,у)/{›* +и* =, умно- 


Уд + В 


жив обе части на —-уАД+В: 


В-уА _ В-уА _ В-уА _ В то —А 
В? -у? А? В?-(1-и’) А? В?-А+и А? В? - А+? А? В-А+и? А?” 
Отсюда 
В —А 
Р= ‚О= 12 
В? — А+ и? А? о В? — А? + и? А? и 


Формулы (12) представляют правило обращения выражения (10). 


Канонические формы ТК-выражений многих переменных можно построить, испо- 
льзуя рекурсивные представления. Пусть Х = {х,,...,х,} < Гама и х, <...<х,. Через 


ЕЁ, обозначим поле ТСоеф ‚ а через Е, (х) — поле рациональных тригонометрических вы- 


ражений от переменной х над полем ГР. Тогда возрастающая последовательность полей 
РА ый СЯ Л 


=Р,(х,.) определяет алгоритм построения канонических форм в поле 


в которой РЁ’. 


рациональных тригонометрических выражений многих переменных. 


3 Определение поля коэффициентов ТВ-выражений 


Выше мы определили поле ТСое/ как поле коэффициентов Т-полиномов. 
В простейшем случае ТСоеГ = Ка. Однако даже в простых примерах тригономе- 
трических задач школьного типа коэффициентами могут быть квадратные радикалы: 
За (л/4) = /2 /2, $ (л/3) = 3/2. 

Использование поля Ка квадратных радикалов уже позволяет расширить класс 
тригонометрических задач, поддерживаемых программной системой. Однако наиболее 
приемлемым является расширение поля Аа тригонометрическими функциями 
от числовых аргументов вида Ал/и,К,пе №. Обозначим через ТСоеГ поле 


Вад (З1щ(К л/п), Соз(Кл/п), К, п Е №). Поскольку 5м(а) = Соз(л/2-@), 
ТСое{ = Кад(Соз(К ли, К, п Е М). (14) 


Поскольку число Со5(Кл/п) алгебраично, для Со5(л/п) над полем Кай су- 
ществует минимальный полином, то есть полином Р(х) е Аа х] минимальной степени 
такой, что Р.(Со5(л/п)) = 0. Преобразованиями понижения степени вида 

Со5" (х) = с. Со5(йх) + с Со5((&-Т)х)+...+ с, Со5(х) + с, (15) 
полиному Р(х) можно поставить в соответствие Т-полином, который естественно 
называть минимальным Т-полиномом числа Со5(л/п). 


Определение 3. Минимальным Т-полиномом числа а = Соз(л /п),п Е № над полем 
Ка! называется Т-полином минимальной степени т вида 
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Т(х) = Соз(тх) + с.Со5((т —Т)х)+...+ с„Со5(х)+с„.с, ЕО, (16) 

такой, что Т(а)=0. 
Следующая теорема устанавливает соотношение между минимальным полиномом 
алгебраического числа е, =с0$(2л / п) +15щ(2л/п) и минимальным Т-полиномом его 


действительной части Ве(е,) = Со5(л/2п). 

Теорема 3. Пусть е, = Со5(2л / п) +195т(2л/п) — первоначальный корень степени 
п из единицы Р,(2) Е Ка{2] — минимальный многочлен е, над полем Кай, Т,(х) — мини- 
мальный тригонометрический полином числа ге(е,) = Со5(2л/п) над Ка и т=4езФ). 


Тогда для любого 2, 


2|=1, имеет место соотношение 


12, 
В (2) =2" °Т,(те(2)). (17) 
Доказательство мы опускаем. 
Следствие. Если Р,(2)=2" +с2"' +...+с„.2+1- минимальный полином перво- 


образного корня е, = Со5(2л / п) +15т(2л/п), т = 21, минимальный Т-полином числа 
Со5(2л/п) имеет вид 


Т(х) = 2(Со5(К) + с Со5((1-Т)х)+...+ с, Со5(х))-+с,. 
Таким образом, для х=2л/п имеет место соотношение 
Со5(2 п/п) =-с.Со5(2(1-Тлп/п-...-с,..С05(27/п)-с,/2. (18) 
Пример. Построение минимального тригонометрического полинома числа 
Со$л/5). Рассмотрим первообразный корень числа е\ = Со5(2л /10) +19т(2л /10). 
Построим минимальный многочлен числа е’. Поскольку это число — корень степени 
10 из 1, его характеристический многочлен равен 2° -1. 
20 —1= (2° +10 (2°-№=(2+0(@“-2+2-2+1(2° -1) 
Отсюда Р‚(х) =х“ -х’ +х? -х+1. Итак, 4ее(Р,)=т=4,[=2. 
То(х) = 2с0$(2х) - 2с05(х) +1 


2с05(2л/5)-2с0$(п/5)+1=0. 


Формула (18) используется для построения канонической формы числового целого 
Т-полинома. 

Определение 4. Числовым тригонометрическим полиномом называется выра- 
жение вида 


Е(Со5(кл /п,), Со5(Кл /п,),..,Со5(К,л/пт„)), Е(жх,...х„)е Ох ,...х„|. (19) 


Определение 5. Канонической формой (стандартным видом) выражения (19) 
называется линейная комбинация вида 


М 
КР)= Ус, соз(}л/п), М <п/2. 
7=0 
Алгоритм вычислений в поле ТСое/ редуцирует числовые Т-полиномы к кано- 


нической форме с помощью (18). Алгоритм построения минимального Т-полинома 
числа со$(27л/п) использует алгоритм построения минимального полинома перво- 


образного корня е, =с0$(2л/ п) +1зш1(2л/п) и (17). 
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Операция деления в поле Тсое/. Выше мы определили алгоритмы операций 
сложения, вычитания и умножения элементов поля ТСоеГ. Мы выяснили, что нетри- 


виальной для реализации является операция умножения двух элементов поля, являю- 
щихся значениями тригонометрических полиномов вида 
Т(х) = Соз(тх) + с, Со5((т-Т)х)+...+с„Со05(х)+с„.с, ЕО 


2л 
при х = —— . Осталось определить алгоритм операции деления двух элементов ТСоеф ‚ то 
Г 


ОВ 
Т,(х) 
ТСоеГ является алгебраическим расширением поля О, мы можем применить стан- 
дартный алгоритм, который состоит в следующем: 

Находим элемент ТСое/ ‚ обратный к Т,(х). Для этого применяем расширенный 


есть алгоритм вычисления коэффициентов полинома Т(х)= . Поскольку поле 


алгоритм Евклида к паре полиномов М (у), Т,(у) , где М(у)- минимальный многочлен 
2 

числа со$(х),х = Е Расширенный алгоритм Евклида вместе с НОД(М(У)Т,(у)) нахо- 
п 


дит такие многочлены О‚(у),0,(у) ‚ которые В} (У)М(у)+1%().0>) = НОД(М(У),Т.(у)). 

Поскольку М(у) неприводимый над О и 4е5(Т(у))<4е(М())), имеем 
НОД(МС), Т,(у)) =1. Итак БО)Ммо)+Ь, ОТ, (>) =1. Поскольку М (х) =0, 
р, (х)Т, (х) =1, или О, (х) = ИТ, (х). 


Этот метод можно модифицировать для тригонометрических многочленов. Для 
этого рассмотрим тригонометрическое тождество 


соз(а)со5(Ь) = соо +Ь) + ооз(а -Б). 


Пусть а=тх,6 =(п-т)х. Тогда соз(тх)со$((и — тух) = сочи) + оз —п)х), 


или 
со5(их) = (2с05(и — т)х))со$(тх) — со5((2т —п)х). (20) 

Равенство (20) используется в расширенном алгоритме Евклида для выравни- 
вания степеней тригонометрических многочленов, старшие члены которых мы обозна- 
чаем через со5(их), со5(тх), т < п. 

Модифицированный расширенный алгоритм Евклида для Т-полиномов. 

Вход: Р, О — тригонометрические многочлены от переменной х. 

Р=сьс05(их) + с, с03((и —Т)х)+...+ с 

О = соз(тх) + с, с0$((т —Т)х) +...+ си 

Выход: О, Т! Т›- тригонометрические многочлены от переменной х. 

Инвариантные соотношения 

р =сСсСБ(Р,О), ТР+Т,О=р 

Система переписывающих правил 

Ех г199Са = г5з(Р, 0,1, \1, 12, \2, М) 


(Р, 0) = (Р,0,1,0,0,1), ИР = 1*Р + 0*О, О =0*Р + 1*@ 
(0,0, 11, \1, 02,\2)=(0, 12, \2), ИО = НОД(Р,О), ТЛ = Ц2, Т2 = \2 
ред(Р) <= 0е9(0) > 

(Р, 0,11, \1, 92, \2)=(0-М*Р, Р, 12-М*И1, \2-М*\1, 11, \1), 
бед(Р) > 0е9(9) > ( 
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(Р, 0,1, \1, 12, \2)=(0, Р, 12, №2, 11, \1) 


7 
Т-моном М (правило 3) вычисляется по формуле 


М = 2% со5((п — т)х), п = 4ее(О), т = дее(Р). 21) 
10 
Отметим еще раз, что алгоритм Ех{Т Сс отличается от расширенного алгоритма 


Евклида для многочленов лишь формулой (21) вычисления дополнительного множителя 
М и интерпретацией умножения как умножения тригонометрических многочленов. 


Выводы 


Методы тригонометрических вычислений, основанные на построении канонических 
форм целых и рациональных тригонометрических выражений и использующие вычис- 
ления в поле коэффициентов ТСоеЁ = Кад(Соз(К л/п, К, п е №) позволяют разрабаты- 


вать эффективные алгоритмы решения стандартных тригонометрических задач — задач 
на упрощения тригонометрических выражений, доказательство тригонометрических 
тождеств и решения тригонометрических уравнений в системах компьютерной матема- 
ТИКИ учебного назначения. 
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М.С. Львов 

Реалзащя тригонометричних обчислень в системах комп’ютерно! математики 

навчального призначення 

У робот! розглянуто алгоритми реалзацй тригонометричних обчислень, основан! на побудовах канон!чних 
форм тригонометричних вираз1в, що застосовуються, зокрема, при програмуванн! стандартних задач 
шюельного курсу тригонометри: задачи спрощення тригонометричного виразу, доведення тригонометрично! 
тотожностт, розв’язання тригонометричного равняння та нер1вност1. 


М.5. Гуоу 

Тиропотейле Сотршайоп$ ш Маетайса! Едисайопа! Зойууаге 

ТБе а2огибил$ оЁ в1еопотеёлс сотрийопз, Базе оп @е сопзйасНоп оЁ сапошса| Юпи$ оЁ @1еопотевлс 
ехргез$1юпз, \ЛисЬ ате ш рагясшаг зе ай ргоотаплплте о запдага 1азК$ оЁ бпеопоглейу зсБоо| сочге, 1.е. 
зпирИйсайопз о ехргезз1оп, ргоуше оЁ14епебу, зоуте оРедианоп ап4 тедиаЩу, аге сопу4егеа п фе ууогК. 


Статья поступила в редакцию 10.06.2011. 
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